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The experimental realization of dynamo excitation as well as theoretical and numerical 
examinations of the induction equation have shown the relevance of boundary conditions 
for a self-sustaining dynamo. Within the interior of a field producing domain geometric 
constraints or varying material properties (e.g. electrical conductivity of the container 
walls or localized high-permeability material) might also play a role. Combining a grid 
based finite volume approach with the boundary element method in a hybrid FV-BEM 
scheme offers the flexibility of a local discretization with a stringent treatment of insulat- 
ing magnetic boundary conditions in almost arbitrary geometries at comparatively low 
costs. Kinematic simulations of dynamo action generated by a well known prescribed 
mean flow demonstrate the reliability of the approach. 

Future examinations are intended to understand the behavior of the VKS-dynamo 
experiment where the field producing flow is driven by ferrous propellers and the in- 
duction effects of conductivity/permeability inhomogeneities might provide the required 
conditions for the measured dynamo characteristics. 

Introduction. Nowadays, there are nearly no doubts that the mechanism 
which is responsible for the generation of astrophysical or planetary magnetic fields 
is a dynamo process in which kinetic energy from a suitable flow of a conducting 
fluid is transfered into magnetic energy. Although the basic idea of this process 
has already been presented at the beginning of the 20th century PQ only few years 
ago fluid flow generated dynamo action has been realized in the laboratory [S1[3J[3]. 
The key parameter that determines the onset of dynamo action is the magnetic 
Reynolds number Rm = /i cVL where a denotes the electrical conductivity, [j,q 
the vacuum permeability, V a typical velocity magnitude and L the characteristic 
size of the considered system. Critical values that are necessary to obtain dynamo 
action in the laboratory are of the order Rm™' ~ 30. ..100 which is already tech- 
nically demanding. Therefore, essential efforts are concentrated on possibilities to 
reduce this critical value and to increase the actual Rm of the field producing flow. 

From numerical simulations it is known that the boundary conditions and also 
boundary layers of stagnant or somehow guided flow could possess supportive as 
well as obstructive impacts on the onset of dynamo action [5, 6J . In the kinematic 
regime the backreaction of the field on the flow by the Lorentz force can be ig- 
nored so that the complexity of the underlying system of equations is significantly 
reduced because only the induction equation with a prescribed velocity field has 
to be solved numerically. Nevertheless, analyzing laboratory experiments requires 
a flexible numerical scheme that is able to consider geometric constraints as well 
as material properties like conductivity jumps between fluid and container walls 
or the high-permeability domains brought in by the iron propellers that are used 
to drive the flow in the VKS experiment, at least in the realization which showed 
dynamo action up to present [3] . The finite volume (FV) approach provides a fast 
and robust scheme relying on a local discretization which delivers an accurate so- 
lution of the kinematic dynamo problem and intrinsically maintains the solenoidal 



character of the magnetic field [7J. However, a drawback of grid based schemes 
are the difficulties arising from non-local boundary conditions as they exist in the 
laboratory in terms of insulating boundaries. Insulating boundary conditions in 
non-spherical geometry in general, are treated by elaborated schemes, e.g. solv- 
ing of the Laplace equation in an extended domain and applying some matching 
conditions [HI [HI QH] , embedding the domain in a sphere [TT] or by simplifying ap- 
proximations (pseudo vacuum, vanishing tangential field). Rather precise results 
which consider insulating boundary conditions exactly are provided by the inte- 
gral equation approach [121 113j . However, the application possibilities are limited 
because of enormous computational resources that are required by this method. A 
different approach that needs less computational power is known as the boundary 
element method (BEM) [14]. Outside the conducting region, the magnetic flux 
density B is expressed as the gradient of a scalar potential B — — V<£, where & is 
determined by the Laplace equation: = 0. Making use of Greens second theo- 
rem and an appropriate discretization, the Laplace equation is integrated only on 
the boundary which in the numerical implementation requires the solution of an 
algebraic set of equations. A combined finite volume/boundary element method 
(FV-BEM) for the induction equation was introduced in [15] . The presented appli- 
cations were restricted, however, to the decay of an initial homogenous magnetic 
field. A deeper investigation of the reliability of the numerical approach in case of 
more complex and realistic problems like dynamo action is still missing. 

The scope of the present paper is the introduction of the methods and basic 
properties of the combined finite volume/boundary element method. The resulting 
scheme is adopted for the numerical solution of the kinematic induction equation in 
cylindrical coordinates in three dimensions with emphasis on the implementation 
of insulating boundary conditions. The integration of further physical effects like 
small scale induction effects parameterized by an a-effect (small scale helical turbu- 
lence) are easily carried out and an extension to a scheme that considers variations 
and/or jumps in conductivity (a) respectively permeability (/i r ) is straightforward 
if corresponding averaging procedures for a or /i r are applied so that the jump 
conditions for E and B at material interfaces are fulfilled [IB]. The detailed de- 
scription of the corresponding methods is beyond the scope of this publication and 
will be presented in a subsequent paper. 



1. Equations and numerical methods. 



1.1. Finite volume method. From Faraday's law dfB = — V x E with the 
magnetic flux density B and the electric field E given by 

1 B 

E = -vxB-\ Vx— (1) 

a flo [ir 

one immediately retrieves the induction equation 

?=Vx (vxB- J-Vx-V (2) 
at V ^Mo Mr/ 

Here, v denotes the velocity field, a the electric conductivity, fiQ the vacuum 
permeability given by fio = 4tt x 10 _7 VsA _1 m _1 and /j, r the relative permeabil- 
ity, that describes the ability of the magnetic flux to penetrate a medium. For 
most substances like air or non-ferrous conducting materials (copper, sodium) 
fi r is very close to 1 whereas ferrous material exhibits a relative permeability in 
the range of ^ r ~ 10 2 ...10 4 . In the following, only homogenous (a = const), 
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non-ferromagnetic (/i r = 1) materials are considered. In a finite volume method 
the computational domain is divided into (small) control volumes where the con- 
servation of all variables is enforced across the control surfaces (interfaces be- 
tween neighboring cells). Writing the induction equation in conservative form 
dtB + V x E = 0, the update of the ^-component of magnetic field at a timestep 
n + 1 in a Cartesian system is given by 
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(3) 



In Eq. ([3]) E* denotes the electric field at an intermediate time step: in an explicit 
scheme with second order accuracy in time * represents the time after the update of 
the predictor step whereas in an implicit scheme * represents the actual timestep 
n + 1. As indicated by the indices, the localization of the components of E is 
slightly staggered with regard to the components of B. Fig. [T] shows the position 
of the field components around a single grid cell labeled (ix, iy, iz). 

Vector quantities are defined on the faces (B, labeled by one half-integer 
index) respectively on the edges of a grid cell (E, labeled by two half-integer 
indices) whereas scalar quantities like conductivity/permeability are defined at 
the center of a grid cell. It is convenient to decompose the electric field into an 
inductive part E md a»xB and a diffusive part E dlS oc (^oO") _1 V x B. E inti is 
treated explicitly applying the C-MUSCL method introduced in [7]. Here, only the 
basics of the scheme are rewritten exemplary for single components of the involved 
quantities. In a predictor step (at an intermediate timestep n + the magnetic 
field on the edges of a grid cell is computed by 
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where the time derivative on the right hand side is computed from the known 
electric field E n at timestep n. The second term on the right hand side dB x /dy is 
approximated using a monotonized central slope limiter which ensures positivity 
preserving, non-oscillating solutions: 



, f Bf y+1 - fl? / Bf -Bf 2 Bf v - Bf y _ 



dy V 2A V ' V Ay Ay ) ) ' 

(6) 

In Eq. ^ minmod(a, b) stands for the minmod limiter defined by 
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< \a\ and ab > (7) 
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The electric field at the intermediate time step n + l /i is then obtained from the 
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upwind solution of a 2D Riemann problem and is given by: 
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where the magnetic field components are the time centered predicted states in- 
terpolated at the edges as determined from Eqs. (j4|) & ([5|). The final update for 
the magnetic field is then performed as described by Eq. ([3]) . In comparison with 
a simple scheme where v and B on the edges of a grid cell are computed apply- 
ing simple arithmetic averages, the C-MUSCL scheme allows for significant higher 
magnetic Reynolds numbers without becoming unstable or exhibiting oscillating 
solutions. 

In an explicit scheme the timestep St is determined by the Courant-Friedrich- 
Lax criteria 

{! = c -^i(^r' <a8 -" )W ) (9) 

with the minimum of the cell extension in x, y or z direction As x _ y ^ z and the 
Courant factor C < 0.5. From expression (0 it is immediately evident that the 
timestep is dominated by the diffusive part oc (As) 2 which in a cylindrical system 
becomes extremely small for grid cells close to the axis. To relax the constraints 
of the time step an implicit solver has been implemented. In a first step an 
intermediate magnetic field B* is computed: 

B* = B cxp - StV x —V x B* (10) 

where £? cxp denotes the magnetic field after the explicit update of the inductive 
part as described above and the diffusive part of the electric field is given in 
discretized form by: 

I i (B y r*. x . If x . . II ■ . /»' •• , . \ 

1 , 1 / ix.iy — -T7 ,iz ix,iy — -^,iz—± ix,iy ,iz — ^ iw,ty—l,iz—-^ J 

cj^o x Ho&\ Az Ay J ' 

(11) 

Similar expressions can be written down for the y- and z-component. 

The resulting algebraic system of equations is solved iteratively for B* using 
a simple Gauss-Seidel method. To remain in the framework of the finite volume 
scheme the electric field at time n + 1 is then computed by 

E diB ' n+1 = — V x B* (12) 
which is used for the final update of the magnetic field according to Eq. ([3]) . 
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1.2. Treatment of the cylinder axis. Although Eq. © describes the field 
update in a Cartesian system, an adoption of the scheme in cylindrical or spherical 
coordinates is straightforward (see e.g. OUS]), essentially by the application of 
the appropriate discretization of the operator Vx = (V''C^ — 9 Z , d z — d r , 1 /rd r (r-) — 
1 /rd v ). However, at r = a coordinate singularity exists that prevents the direct 
computation of BV r=0 , B? r _ x ,^ and B? r _y 2 so that these quantities have to be 
treated in a different way. From the requirement of regularity and uniqueness of 
the solution at r = conditions for the behavior of the magnetic field at the axis 
can be derived. Introducing a decomposition in azimuthal modes 

B = Sft(6 m (r, z, t)e imip ), m = 0, 1, 2, 3, ... (13) 

the behavior of the coefficients b r ^f' z at r = is determined by the following - 
mode dependent - relations: 

F)h z 

m = Q: b^ hi = ^ = 0, 

m=l: V[= f- = |I = o, (14) 
or or 

to>2: b r m = bf H = b z m = 0. 

From these conditions the values of b r ^' z at r = are computed for every mode 
to from the values of the corresponding coefficients close to the axis. The extrapo- 
lation to the axis is based on the radial dependence of b r ^f ,z that follows the most 
general expression for vector quantities close to a cylinder axis given by (see |19j ) 

oo/oo \ oo/oo \ 

B z (r, J) = £ A J2 CLr 2n cos(TO^)+]Tr m £ C£ n r 2n sin(m^) (15) 

m=0 \n=0 / m=0 \n=0 / 

1 oo oo / oo \ 

B r>!fi (r, = C s „r 2 " + £ 1 £ C« „r 2 " cos(to^) 



r 

n— 1 m— 1 \n— 



+ 1 £c£ n r 2 " sm(m^). (16) 



m=l \n— 



where in the numerical realization the polynomial expansion in r is truncated at 
n = 2. 



2. Boundary element method. In case of insulator conditions on the 
boundary the magnetic field is computed by the modified integral equation ap- 
proach presented in [15] . The method depends on the ability to compute the 
normal component of the magnetic field on the boundary by the finite volume 
scheme and a conceptual proximity of the finite volume method and the boundary 
element method (concerning the location of the components of the magnetic field 
on a face centered node). The unknown tangential components of the magnetic 
field at timestep (n + 1) are the result of a matrix operation on a vector composed 
of the normal components of B at the surface of the computational domain. In 
the following a sketch of the scheme is given. 

Insulating domains are characterized by a vanishing current j oc V x B = 
so that B can be expressed as the gradient of a scalar field <P which fulfills the 
Laplace equation: 

B = -V<5 with A<P = 0, <£ -> 0(r~ 2 ) for r -> oo. (17) 
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In a volume fl that is bounded by the surface T Greens second identity for a scalar 
function <P and a test- or weighting function G is written as: 



f f d<P dG 

/ GA<Z> - <PAGdfl = / G— - <P^-dT. (18) 
J J an on 

n r 

If <P = the potential is determined by the integral expression 

*(r) = f G{r,r')^l <P(r') dG{ ^ r>) dT(r>). (19) 
J on On 

r 

where G(r, r') is called Greens function or fundamental solution which fulfills 

AG(r,r') = -6(r-r') (20) 

and is given by 

G(r,r') = (21) 

Furthermore, n represents the direction of the normal unit vector on the surface 
element dT and d /dn is the derivative in the normal direction: d /dn = n ■ V so 
that d n <S> — —B n yields the normal component of B on dr. However, for r € T 
Eq. (JT9J) is not valid since on the boundary <P does not fulfill the Holder criteria 
(|<£(r) — #(t"o)| < Ar a V r with r < c and a, c,A > 0) at every point r, which 
is an essential requirement for (fT5|) . The validity of Eq. (fTTH) can be extended to 
points r G r if the integration domain around a certain source point r located on 
the boundary is enlarged by a small half sphere with the radius e and establishing 
the limit e -> (see Fig. H). 

Writing the enlarged domain r' = V + T e where r e denotes the surface of the 
half-sphere with radius e the first expression on the right side of Eq. (fH?]) becomes: 

lim f G{r, r')^^-dT =lim [gOt, r')^-^dT + lim [G(r,r')^^-dT. 

£^o J on e^o J on a~*o J on 



r'-r. 



With dT e = e 2 cos ddOdtp and \r — r'\ = e the second expression on the right side 
of Eq. (f22|) vanishes: 

lim / G(r,r')^^dr = - lim / -^-^^e 2 costfdddp = 0. (23) 
e^o J dn e^o J 4ire dn 

r e r e 

In the same manner, the second contribution to the integral expression on the 
RHS of Eq. (fl9|) is computed as: 

lim ( dG ^ rl) ${r')dT{r') = lim [ <P(r')dT(r') 
e^o J dn c^o J dn 

r' r'-r,, 

+ lim f dG ^ r '^ (r')dT{r r ). (24) 
c^o J on 
r e 

Unlike in the former case, in Eq. (|24p . the contribution of the integration over 
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the e-sphere does not vanish since 

lim f dG ^ r '^ (r')dr(r') = - lim [ - ' ^ - r 'l $(r')dT(r') 
e^oj dn e^oj 47r|r-r'| 3 

r e 

lim / —^r$(r')e 2 cosfiMdw 

e^O J 47T£ 2 



r 8 r e 



■d = + T/i ip = TT 



J j -^$(r')coadddd<p = -$(r). (25) 



0=— */t <f=0 

Using Eqs. (|22M25[) the integral equation (|19[) is re-written for r £ T as 

^ (r) ./ G(r , r - ) ^2_, M ^ jrM . 

-B n {r') 



(26) 



which is called boundary integral equation. From Eq. (|26p the tangential com- 
ponents of the magnetic field on the boundary B 1 — e T ■ B = —e T ■ V<P(r) are 
computed by: 

B T = 2je T - ($(r')V r dG{ ^ r>) + B«(r')V r G(r, r')) cff(r') (27) 
r 

where e r represents the tangential unit vector on the surface element cff(r'). 
Eq. (|26p and (|2T[) have been derived for a bounded region. An infinite volume is 
treated by introduction of a fictitious surface T describing a sphere with radius R 
in the limit R — > oc. Consider exemplary Eq. (|26p which is re- written including 
terms from infinity: 

i*(r) = | G(r 7 r')^ldT(r') + J G{r,r')^LdT{r') 

^O^W) - / (28) 
an / art 



In three dimensions the following asymptotic behavior of the surface element dT 
and the fundamental solution G(r, r') occurs for R — > oo: 

dT(r') = — 0(7? 2 ), (29a) 

G(r,r') ~ OfE -1 ), ref, (29b) 

~ 0(W\ (29c) 

2 

where J' denotes the Jacobian. Combining the regularity condition oc 0{R ) 
from p7|) at infinity with the asymptotic behavior given in (|29aH29c|) it is ensured 
that the integral expressions in (|28[) that involve R vanish for R — > oo. 

A discretization of the system (|26|) and |27|) yields an algebraic system of 
equations which allows the computation of the (unknown) tangential components 
of the magnetic field. The natural way to define the boundary elements is an 
application of the tessellation provided by the finite volume discretization on the 
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domain surface, where every element has one face centered node at which the 
normal field component B n is located (see Fig. [3j. 

After the subdivision of the surface T in "small" boundary elements Tj with 
T = UTj the potential $i — <P(ri) and the tangential field B\ — B t (ri) = —e r ■ 
(V<Pi) are given by 




Eq. (|30| introduces a global ordering of the quantities <P and B t,n defined by 
an explicit mapping of the grid-cell indices (ix, iy, iz) on a global index (i) with 
i = 0, 1, 2, • • • , N where N — 2 • (nz ■ ny + nz ■ nx + ny ■ nx) represents the total 
number of boundary elements. Then ([>, can be considered as large vectors 
and abbreviating the integral expressions in (f50|) with TL k (k = 1,...,4) the system 
can be re-written in a matrix representation: 

^ = -Hl^j-H^Bf, (31) 

bt = n!'i> s ■ H; s n-. (32) 

Finally, a linear, non-local expression for the tangential field components in terms 
of the normal components results: 

= |ft 3 ® Q ■ 1 - U 1 ^ ®U 2 + H 4 ^j®B n = M®B n . (33) 

The numerical computation of the matrix elements TL^ is performed apply- 
ing a standard 2D-Gauss-Legendre Quadrature method. However, for i — j the 
integral expressions in Eq. ([3H)) become singular so that have to be treated 
separately Since B t is computed from the derivative of <P the scalar potential is 
only determined except an additional constant. Fixing this constant results in a 
relation between diagonal and off-diagonal elements of the matrices Ti. 1 and 7i 3 
(details see [15]): 

The matrix elements H.^ are weakly singular and are computed numerically 
without further difficulties after applying a special cubic coordinate transformation 
where the Jacobian of the transformation has a minimum at the singularity [20| . 

The integral expression that determines the diagonal element of Ti 4 exhibits 
a strong singularity with a vanishing Cauchy principal value so that can be 
computed by excluding some small e- vicinity around r^. 
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Since the elements of 7i k only depend on the geometry and the discretization 
of the problem the computation of M. has to be carried out only once. However, 
since the computation of B t — M. ® B n requires a matrix multiplication with a 
matrix M. of size (27V x N) a large amount of memory is required for M. which 
limits the maximal achievable resolution. This restriction is slightly less severe in 
cylindrical coordinates where the periodicity in azimuthal direction reduces the 
necessary size of the matrix M.. 

3. Results. 

3.1. Simple test case. The free decay of a magnetic field is a simple test 
problem where only the diffusive part of the induction equation (dfB oc AB) 
is considered. A (small) challenge for the scheme arises by initially randomly 
distributed field components which involves the presence of all (resolvable) modes. 
The higher modes decay rather fast and the final solution is dominated by the 
axisymmetric dipole mode which corresponds to the eigenfunction of the system 
with the lowest eigenvalue. Fig. 2] shows the structure of the decaying field in a 
cylinder of height H = 2 and radius R — 1 after approximately one diffusion time 
(Tdifi = MoC-R 2 ) which is dominated by an axial dipole. 

The temporal behavior of the total magnetic energy i? mag = (2/io) -1 / B 2 dV 
is shown in Fig. [5l After the system settles down to its eigensolution (at t ~ 
O.lTdiff) a simple exponential decay oc e - At /^dif f i s 

observed. 

The influence of the boundary conditions is obvious in comparison with the 
decay rate in case of vanishing tangential field boundary conditions (VTF, dashed 
curve) where a significant slower decay takes place. Table Q] shows the decay rate 
for both types of boundary conditions in comparison with the results obtained 
by [H| applying the integral equation approach (IEA) and a differential equation 
approach (DEA) where the Laplace equation is solved in the exterior. 

The decay rate of the m = mode in case of insulating boundaries is larger 
than for VTF conditions because of smaller or even vanishing field gradients in 
the latter case. The decay rate for the m = 1 mode achieved from the FV-BEM 
scheme is approximately 4% larger than the decay rate obtained from the IEA. 
The observed deviation between the different schemes might arise through a couple 
of reasons which will be discussed in the following subsection. 

3.2. Von-Karman Row. The experimental realization of a von-Karman-likc 
flow in the VKS experiment is driven by two counter-rotating propellers located 
close to the end-plates of the cylindrical vessel. An analytic prescribed velocity 
field that roughly resembles the mean flow in the VKS experiment is given by the 
so called MND-flow (Marie-Normand-Daviaud, [21]): 

v r = -0.57rcos(7rz) r(l - r) 2 (l + 2r) 

v v = 4er(l — r) sin (0.5ttz) (35) 

v z = (1 — r)(l + r — 5r 2 ) sin (nz) 

where e describes the relation between toroidal and poloidal component of the 
velocity (here given by e — 0.7259 which has turned out to be an optimum value 
for generating a dynamo |21|). 

The structure of this flow is essentially dominated by two toroidal and two 
poloidal cells which are visualized in Fig. [S] where the grey scaled contours rep- 
resent the azimuthal flow component and the arrows represent the poloidal flow 
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components. If the Reynolds number - defined as Rm = a^ioRV m&x - exceeds a 
critical value, dynamo action takes place. Detailed examinations of the induction 
effects of the MND-flow - in particular on the effects of side and lid layers contain- 
ing a stagnant fluid component that surrounds the active region - can be found in 
[BJ. Here, only simulations without outer fluid layers are discussed. 

The initial condition is again given by a divergence-free random field, which 
- if the applied magnetic Reynolds number exceeds the critical value - after a 
timespan of about a diffusion time passes into a growing smooth field (Fig. [7]) 
field. The obtained field structure essentially represents the well known banana- 
cell like behavior (|6j|22]) corresponding to a field dominated by the m = 1 mode 
(see left hand side of Fig. [7j). The right hand side of Fig. [7] presents a snapshot 
of the streamlines of the magnetic field from a typical simulation (Rm = 80) from 
which the clear dominance of the equatorial oriented dipole solution is apparent. 
A quantitative comparison between the FV-BEM and the IEA schemes has been 
done by means of the growth-, respectively decay rates for the m — 1 mode. 

If Rm Rm crl the field decays and the slowest decaying eigenmode can be 
observed. For the 3D grid based FV-BEM scheme the growth rates have been 
computed from the temporal behavior of the volume averaged Fourier amplitudes 
(azimuthal decomposition) of the axial field componenfl The growth rates of the 
dynamo (or decaying) state in dependence of the applied Rm are presented in 
Fig. U 

Both schemes provide rather similar results for small and moderate Reynolds 
numbers. Small but systematic deviations between the IEA and the FV-BEM 
occur for Rm ^ 45. The main and most probable reason for this characteristic 
behavior might be the rather low resolution attainable by the IEA which is re- 
stricted to 20 x 20 points in a 2D simulation due to the enormous computational 
power required by this method. A further influence might be obtained from the 
coupling of different field modes that necessarily occurs in the 3D FV-BEM scheme 
whereas the IEA only considers one single mode. However, it seems unlikely that 
this influence has a remarkable effect since the amplitudes of the non-dominating 
modes are always suppressed by several orders of magnitude. 

For the FV-BEM scheme the critical Reynolds number is obtained from in- 
terpolation of the growth rates around the occurrence of dynamo action and is 
given by Rm crlt = 62.6. This value is situated slightly above the value reported in 
[22] (Rm crit = 58) or the results of simulations with the DEA (Rm crit = 61.5) re- 
spectively the IEA (Rm crlt = 59.6) reported in |6j. For idealizing VTF conditions 
the critical Reynolds number is computed as Rm clIt = 41.2 (see dotted curve in 
Fig. [8] The significant deviation between the critical Reynolds number for phys- 
ical boundaries and for idealizing (VTF) boundaries depends on the geometry of 
the system and becomes smaller for a larger relation between height and radius. 
Beside Rm c " , a second distinguished point exists, at which the transition of a 
m = dominated solution to a m = 1 dominated solution proceeds, which occurs 
at Rm « 25. Below this value the axial dipole dominates the field structure, how- 
ever, the m = mode never becomes unstable to dynamo action. Note, that for 
small Reynolds numbers both schemes yield a slight minimum for the decay rate 
of the ni — 1 mode, so that for Rm ^ 15 this mode decays faster than without 
any flow. 

4. Conclusion. A fast and easy to handle code base for 3D simulations 
of the kinematic induction equation has been developed which is able to treat in- 

1 All components of the magnetic field behave similar, however, B z is the dominant component. 
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sulating boundary conditions in a stringent way. The reliability of the approach 
has been tested against known results of dynamo action generated by the (analytic 
prescribed) MND-flow and a good agreement was achieved in comparison with the 
results of the integral equation approach. The IEA is also able to treat insulating 
boundary conditions but requires much more computational power and provides 
less flexibility concerning the addition of further physical effects or geometric con- 
straints. 

Significant differences for the critical Reynolds number and the growth rates 
exist between simplifying boundary conditions and insulating boundaries. These 
discrepancies become smaller for higher aspect ratios. However, in case of a real- 
istic relation between height and radius as it usually is realized in the laboratory 
a thorough consideration of the appropriate boundary conditions is indispensable. 

The flexibility of the scheme facilitates the addition of further physical terms 
like an a-effect as additional dynamo source or the consideration of temporal 
fluctuations of the mean flow. A direct extension of the scheme which will be 
presented in a subsequent paper will be able to consider conductivity/permeability 
inhomogeneities. 
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Figure 1: Localization of the components of the electric/magnetic field on the 
Cartesian staggered mesh. Not all components are labeled for the reason of clarity. 




Figure 2: Enlarged domain around the source point r located on the boundary 
surface T. The boundary integral equation is evaluated over an extended domain 
r + T e in the limit e -> 0. 
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Figure 3: Grid structure and boundary element discretization. i,j denote a global 
ordering number. 
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Figure 4: Free decay of an initially random magnetic field. Snapshot of the 
field structure after one diffusion time. Note the non-vanishing contributions of 
the tangential components at the boundary caused by the insulating boundary 
conditions. 
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Figure 5: Growth rate of a freely decaying magnetic field in a cylinder. The solid 
(dashed) curve shows the total energy applying insulating (vanishing tangential 
field) boundary conditions. 
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Figure 6: Toroidal (shaded contours) and poloidal (arrows) components of the 
von-Karman-like flow prescribed by Eq. (135I) . 



16 





Figure 7: Left panel: Isosurface of the magnetic energy at 20% of the maximum 
value. Right panel: time snapshot of the magnetic field structure represented by 
the streamlines. Rm = 80. 
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Figure 8: Growth rate of the dominating modes (m = 0, grey solid curve and 
tyl — 1 j black solid curve) in dependence of the magnetic Reynolds number Rm. 
The black dashed line corresponds to the growth rates of the m = 1 mode obtained 
from the IEA [6] . The dotted curve shows the results for vanishing tangential field 
conditions (VTF). 
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